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Introduction. In ultrarelativistic heavy-ion collisions, the pion freeze-out phase-space 
'k>( I density determines the importance of multiparticle pion correlations and of dileption pro- 
;h ' duction from tt+tt" annihilation. G. Bertsch suggested a way of estimating this quantity 
by combining measurements of one-particle momentum spectra with two-particle correla- 
tions, thereby testing the thermal equilibrium of the pion source created in the collision. 
In local thermal equilibrium at temperature T{x) the pion energy distribution is given by 
the Bose-Einstein function 

J[X,P) = ^p.u(x)/T(x) _ I ' ^1) 

where u{x) is the 4- velocity of the local rest frame at point x in the observer frame. The 
coordinate space average of this function is the quantity to be measured: 



J f{x,p)pf'd^af,{x) 



Here d^a{x) is the normal vector on a space-like space-time liypersurface a{x). According 
to Liouville's theorem, a is arbitrary as long as its time arguments are later than the time 
if(x) at which the last pion passing the surface at point x was produced. 

If the measured single-particle p^-spectrum is parameterized by an exponential with 
inverse slope parameter Tef[{y), 

dn- _ dn- 1 / Pt \ ,^. 

dyprdpT dy Tjj(y) ^ %a{y)) ' 

and the two-particle correlation function by a Gaussian 0, 

C{q] Pt, y) = 1 + A(pt, y) exp - ^ RI{pt, y)ql - '^Rliipr, y)qoqi , (4) 

y i=o,s,l J 

where the subscripts o, s, I refer to the usual out-side-long coordinate system and \{pT-, v) 
accounts for unresolvable contributions from long-lived resonances , one obtains for the 
spatially averaged phase-space density at freeze-out [§|, |1 

{f){PT,y) = JM^) '' 'r^^^'^ (5) 

7r-3/2 Ep R,{pT, y)jRl{PT, y)RKPT, y) - RtiiPr, y) 



Here Ep=\/rnF+^=mT coshy, with rriT = Jrn? + p"^. The numerator (with experimental 
input dn~ /dy,Teff{y)) gives the momentum-space density at freeze-out while the denom- 
inator (involving the measured two-pion Bose-Einstein correlation radii Rs,Ro,Ri,Roi) 
reflects the space-time structure of the source at freeze-out and can be interpreted as 
its covariant homogeneity volume for particles of momentum p. The factor vA ensures 
1^ that only the contributions of pions from the decays of short-lived resonances, which 
happen close to the primary production points, are included in the average phase-space 
density. 

We have calculated {f){pT,y) for the S-S, S-Cu, S-Ag, S-Au, S-Pb and Pb-Pb data 
from the experiments NA35 [^ , NA49 []^, ^ , and NA44 |Q , and for the vr-p data from the 
NA22 experiment |]10| (all at the CERN SPS). The projectile energies were 200 GeV per 



nucleon in S-nucleus collisions, 158 GeV per nucleon in Pb-Pb collisions and 250 GeV in 
TT-p collisions, which correspond to projectile rapidities 6, 5.8 and 8.2, respectively. Our 
results will also be compared with the average phase-space density in Au-Au collisions at 
projectile momentum 10.8 GeV/c, published by the E877 collaboration at the AGS 0. 
In all cases the analysis was done in the LCMS (longitudinally comoving system) where 
the longitudinal momentum of the pion pair vanishes. 



Calculation of the phase-space density. The experimental input into Eq. (^) is partly 
incomplete. This concerns particularly the intercept parameter A which was measured 
in all the considered experiments but not always published in the final corrected form. 
However, the necessary information is still available, and in the following we explain how 
we used it. The following experimental effects are responsible for the uncertainties in the 
measurement of the A parameter: 

(i) Finite momentum resolution reduces the correlation intercept, since it leads to a 
smearing or widening of the correlation peak. 

(ii) Corrections for Coulomb repulsion of like-sign pions play an important role in the 
Bose-Einstein correlation analysis, in particular in the measurement of A. Certain 
sophisticated Coulomb correction methods have been applied p, H, |^, [11], |12[, but 
not to all the data sets, as we shall discuss below. 

(iii) If pions are not positively identified, as is the case in NA35 and NA49, the direct- 
pion sample is contaminated with kaons, converted electrons, protons and other 
particles. NA35 and NA49 have performed a contamination correction based on de- 
tailed Monte Carlo simulations [H, but unfortunately the resulting corrected results 
for A have not been directly published. However, since the contamination levels 
were published, as well as the uncorrected results for A, we were able to estimate 
the corrected values ourselves. For example, in the NA49 Pb-Pb collisions, in the 
rapidity interval 3.4< y <3.9, the fraction of pure 7i~7i~ pairs is a: ~ 55% 0, and the 
uncorrected measured Xmeas is between 0.4 and 0.5 0. The contamination-corrected 

value is then Xcorr = Xmeas/X ~ 0.73-0.91. 

(iv) Pions originating from weak decays also contaminate the direct pion sample when 
the decay takes place unresolvably close to the main interaction point. In the rejec- 
tion of decayed pions, the experiments NA35 and NA49, with continuous tracking 
detectors covering the region close to the target, are in an advantage over the NA44 
spectrometer with tracking devices placed only far away from the target. Since 
there are more hyperons (A, S, S) than antihyperons in the fireball which produce 
more negative than positive pions, the measurement of A with positive pions is less 
affected and should be more reliable. Note that NA44 indeed reported P] signifi- 
cantly different results: A ~ 0.52 ± 0.03 and A^+ ~ 0.66 ± 0.04, the difference 
being mostly reproduced in a detector simulation with the RQMD event generator, 
and thereby traced to the weak decay asymmetry. 



NA44 and NA49 both reported an approximate independence of A on p^, valid to about 
±10%. Regarding the absolute values, for S-Pb data NA44 has reported A"'"+(NA44) ~ 
0.52 — 0.59 (with negligible non-pion contamination), consistent with the lower half of 
the NA35 interval, A (NA35) ~ 0.55 — 0.7 (contamination corrected). The latter may 
be read out directly from the correlation functions presented in Fig. 2 of Ref. 0. That 
figure also provides an insight into the reliability of the A measurement, since it presents 
the evolution of the correlation function through different correction steps. 

In case of Pb-Pb collisions the NA44 results are again consistent with the lower end 
of the NA49 interval, see Fig. [l|: A (NA49)~0.7-0.9, obtained with "our" correction for 
contamination, and A^+(NA44)~0.66-0.69, without a correction for the contamination 
due to weak decays. Since positive pions are also contaminated by weak decay products, 
the NA44 result is certainly an underestimate, unlike the NA49 result which has been 
corrected for all sources of contamination, including those weak decays which could not 
be filtered out by vertex cuts. 

If only the contamination by weak decay products were to explain the difference be- 
tween A (NA44) and A"'""'"(NA44), as suggested by NA44 [^, then the contamination 
effect itself on A"'""'"(NA44) should be quite strong, probably on the order of 10%. After 
shifting the NA44 results up by this ad hoc 10% correction, as shown in Fig. |I| by full 
lines, the NA44 and NA35/NA49 intervals already fully overlap. 

In comparing S-nucleus and Pb-Pb results one has to take into account also a difference 
in the Coulomb correction methods. Sophisticated Coulomb correction methods lead to 
higher A measurements than the "old" Gamov correction. This effect seems to be around 
5-15% for Pb-Pb data P, |], and should be lower for S-nucleus data |l^, |ll|, probably 



0-10%, depending on the way the normalization has been done 0, 11, 12|. Note that the 
Pb-Pb data have been Coulomb-corrected in the appropriate way, while the S-nucleus data 
were only Gamov corrected, which might be the reason for the discrepancy seen in Fig. |l|. 
Assuming that the S-nucleus data indeed need to be corrected by 10% to account for the 
systematical error due to the inappropriate Gamov correction, one arrives at the results 
presented by dotted lines in Fig. |I| which are in much better agreement with the Pb-Pb 
results. Our assumption is, however, rather doubtful and should serve as an illustration 
of the systematic uncertainties rather than as a quantitative result. Note that also the 
correlation radii should in principle increase due to a proper Coulomb correction, which 
would partly cancel or even reverse the effect of an increased A, see Eq. (||). 

It is beyond our ability to check all these points in further detail, since this would 
require detailed Monte Carlo simulations of the detector response, with a realistic event 
generator as input. It is quite obvious that with the quality of the presently available data 
an accurate determination of the factor vA to be used in Eq. (^ is not possible. We have 



therefore proceeded in a rather pragmatic way, taking the most probable interval for A 
for all the data sets from Fig. |l|, based on a consideration of the corrected data presented 
by full (and to some extent also dotted) lines in Fig. |I]. The most likely value for all these 
data sets seems to be A= 0.7, with a relatively wide error ±0.1. This value has been used 
in our analysis. 

Fortunately, A enters Eq. (^ under a square root, and a systematical error of e.g. 15% 
for A will result in only 7% error of the final result, which raises the confidence into our 
simplistic approach. 

In the introduction we stated already that we shall use an exponential function, Eq. @, 
to parametrize the pT-spectrum. One of the reasons for our choice was the general ac- 
cessibility of Teff for all the data sets; Tes equals {pt)/^, a quantity which is usually 
quoted together with the data. The best method would be to use directly the measured 
dn~ /(PpX) but most of the data have not been published in a useful form. The pion rapid- 
ity density dn~ /dy was derived from the NA35 and NA49 measurements of the negative 
hadron rapidity density dn^~ /dy, by a simple scaling with a factor 0.9 ||13|| . 



One should also specify the px value representative for a given p-r-interval, both in 
order to calculate (/) and to allow for a meaningful model comparison, i.e. to properly 
place the calculated points in the plot (/)(pt)- Fortunately, the result is almost insensitive 
to systematical errors in the choice of the representative px value. To estimate the effect 
of this error we have shifted each measured point by 10% up and down in pt from the 
average px for the interval, as shown in Fig. ^ The resulting smearing is elongated along 
the actual shape of the entire distribution. The systematical error in the choice of the pt 
value representative for a given p^ interval is therefore essentially harmless. 

With the S-nucleus collision data from NA35 and NA44 and with the vr-p collision data 
from NA22 we had yet another problem: the two-pion correlation functions of these data 
have not been fitted with the proper functional form which includes the out-longitudinal 
cross term Roi [0. For those data the resulting volume term in the denominator of Eq. (^ 
reduces to RsRqRi- To estimate the systematic error arising from the omission of the cross 
term we considered the ratio of the expressions for (/) with and without the cross term: 

<^^ - ' (6) 



{f){Roi = ^) U _ K, 



With the explicit expressions given in f^ and the Cauchy inequality one shows that 
R'^i < RgRf such that this ratio is always well-defined. It must be 1 at midrapidity since 
there Roi = by symmetry. Larger values of 1.2-1.3 were only found for data far away 
from mid-rapidity for which the measured cross term was comparable to Ro or Ri. In this 



publication we have considered data for which no cross term was included in the analysis 
only near mid-rapidity where Roi is small. We estimate that that this procedure limits 
the (downward) systematic error on (/) due to this effect to less than 5%. 

Results. In Figures ^^ we show the average phase-space density as a function oipT, for 
different systems and rapidity bins, as extracted from Eq. (^. The error bars reflect the 
statistical errors of the single-particle spectra and correlation radii, but not the systematic 
uncertainties discussed above (including the dominant uncertainty in the intercept A) . We 
also plot for comparison the Bose-Einstein distribution (|l|) for a static system [u{x) = 0) 



and different temperatures T. In doing so we set E ^ Jm? + p^ which is justified if the 
longitudinal pair momentum can be neglected. Since the correlation data which we used 
were analyzed in the so-called "fixed longitudinal co-moving system" (FLCMS) in which 
the pair rapidity corresponding to the center of the rapidity bin vanishes, this approx- 
imation holds as long as the considered rapidity interval is narrow. Rapidity intervals 
of 0.5 (as in the NA49 Pb-Pb data) and 1 (as for the NA35 S-nucleus data) result in a 
systematic overestimate of the theoretical distribution function by approximately 2% and 
7%, respectively. 

From the results presented in Figures ^|-^, one may draw the following conclusions: 

(i) Universal average phase-space density at freeze-out: 

Even though the heavy-ion data span about an order of magnitude in multiplic- 
ity density (c/n-/rfy(S-S)=22, rfn-/rfy(S-Cu)=24, rfn'/c/y(S-Ag)=28, dn-/dy{S- 
Au)=37, and (in"/(i|/(Pb-Pb)=40-185, depending on rapidity), the resulting average 
phase-space densities vary by much less. Given the large error bars, all the nuclear 
collision data at mid-rapidity from the SPS experiments in Figs. |^,| are almost in- 
distinguishable. Freeze-out happens in all cases at comparable values of the average 
phase-space density. 

(ii) Rough agreement with the Bose-Einstein distribution: 

The dashed lines in Fig. |] indicate thermal Bose-Einstein distributions for static 
sources at various freeze-out temperatures. A rough comparison with the SPS heavy- 
ion data indicates consistency with temperatures in the range of 100-140 MeV. 
Thermal freeze-out temperatures in this domain were recently obtained in analyses 
of the measured spectra and correlation functions H, [11, O |16[. 



(iii) Multiparticle symmetrization effects are small: 

Irrespective of the collision system and pt, we find phase-space densities smaller than 



0.5. For (/) <^ 1 the Bose-Einstein phase-space enhancement is dominated by two- 
particle symmetrization effects, and muhiparticle symmetrization effects are weak 
1^. This is an important consistency checlc for the current practice of calculating the 
two-particle correlation function from the two-particle symmetrized contributions 
only. For (/) < 0.5 the system is far away from a pion laser. We find no sign for a 
striking pion excess or a pion condensate. 

(iv) Radial flow: 

Looking in more detail, in particular applying a logarithmic scale as in Fig. |, one 
finds that the data indicate a somewhat slower decrease with increasing pt than the 
Bose-Einstein curve. Particularly strong differences are seen for the NA49 Pb-Pb 
data at mid-rapidity (2.9 < y < 3.4). These data can be fitted with the function 
exp[-l.l(l) -pT/0.31(6)GeV] (xVndf=0.2), presented in Fig. | by a dotted line, 
which has a considerably flatter slope than the Bose-Einstein distribution. This 
behaviour can be reproduced by a model for the emission function which includes 
radial collective expansion and whose parameters are adjusted to reproduce the 
single- and two-particle spectra |[16|. The (strong) radial expansion adds extra 
transverse momentum to the particles (Doppler blue shift), i.e. the local (/) values 
appear in the lab frame at a higher pt than in the rest frame of the effective source. 



A detailed study will be published elsewhere jlB 



Decoupling at high temperature in 7r-p collisions: 

In contrast to freeze-out in nuclear collisions which takes place in two steps (chemical 
freeze-out of particle abundances at around Tchem — 170 — 180 MeV [|l7l, thermal 
freeze-out of momentum spectra at around Ttherm — 120 MeV 0), pion production 



in vr-p collisions |jTO[ is essentially immediate, without the second evolution stage, 
and therefore common chemical and thermal freeze-out temperatures of around 
170-180 MeV should be expected. The data on (/) are indeed consistent with this 
expectation, as seen in Fig. ^ 

(vi) Rapidity dependence: 

A certain departure from the universal scaling is seen for the data at rapidities close 
to the projectile rapidity, see Fig. |^, both at AGS and SPS; but again the results 
at these two widely different beam energies are mutually consistent and agree with 
the expectations from a thermalized expanding source. The sources at AGS and 
SPS energies show strong longitudinal expansion at freeze-out, but with maximal 
longitudinal flow rapidities well below the beam rapidity (r/flow.max ~ 1-7 at the SPS 
and ~ 1.1 at the AGS |T8|, |19[). Pions with CMS rapidities larger than r/Aow.max 



thus come from the tail of the thermal longitudinal momentum distribution [Q . A 
similar decrease of (/) which is seen near midrapidity as a function of pt is thus 
seen near beam rapidities as a function of y. 
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Figure 1: Correlation strength (intercept parameter) A. Dashed hnes: raw uncorrected 
results. Full lines: contamination corrected results. Dotted lines: ad hoc 10% correction 
for inappropriate Gamov Coulomb correction. 
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Figure 2: Each data point was placed at three different places in p^'- at the average 
value pT of the corresponding bin (sohd), at 0.9 pt (open), and at 1.1 Pt (grey symbols). 
The bin positions are indicated by larger tics. The resulting smearing follows the general 
shape of the distribution. 
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Figure 3: Phase-space density as a function of pt- Different heavy- ion data sets from 
the SPS are indistinguishably similar. Dashed hues indicate Bose-Einstein distributions 
for three choices of the local freeze-out temperature: 80 MeV, 120 MeV and 180 MeV. 
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Figure 4: Logarithmic representation of the data shown in Fig.|^. Discrepancies between 
the data and the Bose-Einstein distribution (T =120 MeV) can be explained by radial 
transverse flow. The dotted line represents an exponential fit to the mid-rapidity Pb-Pb 
data (see text). 
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Figure 5: NA22 7r-p data. They are consistent with a Bose-Einstein distribution of 
temperature T = 180 MeV. 
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Figure 6: Average freeze-out phase-space density for A-l-A collisions {A ^ 200) as a 
function of pt in different rapidity intervals. At very forward rapidities (/) decreases 
both at the SPS and the AGS as expected from a thermalized, longitudinally expanding 
system (see text). 
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